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and four-wave mixing in an optically dense atomic vapor 
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We experimentally and theoretically analyze the propagation of weak signal field pulses under 
the conditions of electromagnetically induced transparency (EIT) in hot Rb vapor, and study the 
effects of resonant four- wave mixing (FWM). In particular, we demonstrate that in a double- A 
system, formed by the strong control field with the weak resonant signal and a far-detuned Stokes 
field, both continuous-wave spectra and pulse propagation dynamics for the signal field depend 
strongly on the amplitude of the seeded Stokes field, and the effect is enhanced in optically dense 
atomic medium. We also show that the theory describing the coupled propagation of the signal and 
Stokes fields is in good agreement with the experimental observations. 

PACS numbers: 42.50.Gy, 32.70.Jz, 42.50.Md 



I. INTRODUCTION 

In recent years, electromagnetically induced trans- 
parency (EIT) [H |2] and the associated effect of ultra- 
slow pulse propagation ("slow light") [3] have attracted 
considerable attention due to their many promising ap- 
plications. For example, they enable coherent, reversible 
transfer between the quantum states of an electromag- 
netic field and the collective excitation of an ensemble of 
long-lived radiators (e.g., spins of atoms, solid-state im- 
purities, quantum nanostructures, etc.), which are nec- 
essary for the realization of quantum memory [3] , entan- 
glement of distant matter nodes [5j , single photon gener- 
ation [6] , and realization of deterministic two-qubit gates 
for photons [7J, etc. At the same time, more robust con- 
trol of the propagation of classical pulses using EIT has 
been extensively explored in various materials for appli- 
cations such as optical packet switching and optical signal 
processing [5]. 

In a traditional EIT scheme, a strong classical con- 
trol field is applied to one optical transition, resulting in 
a modification of the optical properties of a weak signal 
field, which couples the same excited electronic state with 
a second long-lived ground-state sublevel, thereby form- 
ing a A system pQ, as depicted in Fig. [lja). In the limit 
of low optical depth, it is sufficient to take into account 
only the effects of this single A. However, many appli- 
cations require operation at high optical depth [51 HO] , 
where additional nonlinear effects may become impor- 
tant [TTJ [T2 [T3J E]. One such effect is resonant four- 
wave mixing — a nonlinear process arising from the far 
off-resonant interaction of the control field. Earlier stud- 
ies [T5J [TB] found that the propagation of the signal field 
in this case will be strongly affected by the presence of 
the Stokes field. 

The effect of four-wave mixing can be advantageous 
or detrimental, depending on the details of the appli- 
cation. For example, in quantum memory applications 
the resonant mixing reduces the fidelity by adding extra 
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FIG. 1: (Color Online) (a) The double-A system used in the- 
oretical calculations, (b) A schematic of the experimental 
apparatus (see text for abbreviations). 



noise into the signal field. Also, FWM may limit the 
storage efficiency at higher optical depth [TTJ. On the 
other hand, non-classical correlations between two signal 
and Stokes fields can individually carry quantum infor- 
mation |18j and produce entangled images. Similarly, for 
slow light applications, the conversion of an original pulse 
from the signal to Stokes channel may reduce the read- 
out efficiency |19j . However, under certain conditions, 
FWM may lead to gain for both the signal and Stokes 
fields, which could compensate for any optical losses 

In this manuscript, we investigate the modification of 
the two-photon resonant transmission peak and the dy- 
namics of pulse propagation in the case of a seeded Stokes 
field, which was produced simultaneously with the signal 
field by phase modulating the original monochromatic 
control field. Specifically, we show that the resulting sig- 
nal and Stokes spectra are well-described by a double-A 
system. Additionally, we show that this model accurately 
portrays the dynamics of signal and Stokes pulse propa- 
gation. 



II. EXPERIMENTAL ARRANGEMENTS 

Measurements were performed using the configuration 
i Fig. [TJb). We tuned a commercial external cavity 
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diode laser (ECDL) near the Rubidium D\ transition 
(A = 795 nm). After using a polarizing beam splitter 
(PBS) to separate a fraction of the light for a reference 
beam, we passed the main beam through an electro- 
optical modulator (EOM) , which modulated its phase at 
the frequency of the ground state hyperfine splitting of 
87 Rb (A hf /27r = 6.835 GHz). Due to the small driv- 
ing amplitude, the phase modulation mostly produced 
two first modulation sidebands at ±Ay of nearly equal 
amplitudes and opposite phases. We tuned the zeroth 
order (carrier frequency) field to the 5 2 S!/ 2 F = 2 — > 
5 2 P!/ 2 F' = 2 transition; this beam acted as the con- 
trol field. The +1 modulation sideband of the amplitude 
functioned as the signal field, and was tuned near the 
5 2 S 1/2 F = 1 -> 5 2 P 1/2 F = 2 transition. The -1 side- 
band acted as the far-detuned Stokes field. Then all op- 
tical fields passed through an acousto-optical modulator 
(AOM) operating at 80 MHz, which shifted the frequen- 
cies of the fields by that amount. 

For spectral measurements, the control field was al- 
ways on, thereby ensuring that most of the atoms were 
in |g), and we slowly swept modulation frequency of the 
EOM, which synchronously scanned the two-photon de- 
tuning for the Stokes and signal fields. For slow light 
measurements, we first applied a pulse of the control field 
to optically pump the atoms into |g), and then we ad- 
justed the modulation power of the EOM and AOM to 
produce desired Gaussian pulses of the signal and Stokes 
fields. 

To carefully evaluate the effects of resonant four-wave 
mixing, we used a temperature tunable Fabry-Perot 
etalon (FSR = 20 GHz, finesse w 100) to reduce the 
Stokes field amplitude by tuning the etalon transmission 
resonance such that most of the — 1 sideband is transmit- 
ted, but all the other fields are reflected. This way, we 
were able to reduce the intensity of the Stokes field by a 
factor of 20 without noticeable losses in the intensity of 
both the control and signal fields. 

Before entering a vapor cell, the beam was weakly fo- 
cused to either 2.6 mm or 4 mm diameter, as we indicate 
below, and circularly polarized with a quarter-wave plate 
(A/4). Typical peak control field and signal field powers 
were approximately 19 mW and 50 /xW, respectively. A 
cylindrical Pyrex cell, of length 75 mm and diameter 22 
mm, contained isotopically enriched 87 Rb and 30 Torr Ne 
buffer gas, so that the pressure-broadened optical tran- 
sition lincwidth was 2-y = 2tt x 290 MHz [2U- The cell 
was mounted inside tri-layer magnetic shielding, as to 
reduce the effects of stray magnetic fields. The tempera- 
ture of the cell (and thus the concentration of Rb in the 
vapor phase) was adjusted using a bifilar resistive heater 
wound around the innermost magnetic shielding layer in 
the range between 70° C and 80° C, which corresponded 
to the change in Rb densities from 5.6 x 10 11 cm~ 3 and 
1.2 x 10 12 cm~ 3 , and to the range of optical depths 2d 
between 52 and 110. Here we define the optical depth 
2d such that the probe intensity without EIT is attenu- 
ated by e~ 2d , and our procedure for calculating the ef- 



fective optical depth is described in Ref. [IT] . We also 
measured |17j the typical spin wave decay time to be 
l/(27 Sfl ) ~ 500 /is, most likely arising from small, un- 
compensated remnant magnetic fields. 

After the cell, the output laser fields were recombined 
with the unshifted reference beam at a fast photodetec- 
tor, and the beat note signals between each of the +1 
and — 1 modulation sidebands and the reference field was 
measured using a microwave spectrum analyzer. Because 
of the 80 MHz frequency shift introduced by the AOM, 
the different beat note frequencies of each sideband with 
the reference field allowed for independent measurement 
of their amplitudes. 



III. THEORETICAL MODEL 

In a traditional three-level system, under the con- 
ditions of electromagnetically induced transparency, a 
strong control field (frequency vq = u es , Rabi frequency 
f2) |22j works in conjunction with a weak signal field {y = 
LOgg+S, Rabi frequency a) to create a long-lived coherence 
between states |g) and |s) [J, as depicted in Fig. [TJa), 
producing a narrow symmetric transparency peak in the 
spectrum of the signal field near the two-photon reso- 
nance (u — Vq ~ Ugs)- Simultaneously, the signal field 
experiences a steep variation in refractive index, thereby 
reducing the pulse's group velocity ("slow light"), and 
leading to a pulse delay time of r = (d"f)/\£l\ 2 3> L/c [3J, 
where 7 is the optical polarization decay rate. 

While this description is sufficient in the limit of low 
optical depth |17j . at high optical depth, accompany- 
ing nonlinear processes become important. A sufficiently 
strong control field can excite the far off-resontant transi- 
tion from |g) to |e) [15], and generate a Stokes field a' via 
the four- wave mixing (FWM) process. We theoretically 
model this effect by considering a double-A configuration, 
consisting of a near-resonant A system formed by the con- 
trol and signal fields and of the additional far-detunied A 
system formed by the same control field (i>c — ^es, Rabi 
frequency fl'), applied to the state |g) and by an addi- 
tional Stokes field (u 1 — uj es — Ahf — 5, Rabi frequency 
a'). In such system, we can use Floquet theory [33J to 
adiabatically eliminate the off-resonant interaction via f2' 
and a' . As a result, to first order in 1/Ahf and in a', one 
obtains an effective Rabi frequency Wa'* /Ahf coupling 
|g) and |s), while the states |e) and |g) acquire small 
light shifts S s = |f2'| 2 /Ahf and — S s , respectively. For our 
Clebsch-Gordan coefficients [T7], |fi'| 2 = 3|fi| 2 . 

The rotating-frame Hamiltonian describing this inter- 
action is: 



H = -(S-S s )\s)(s\-(S-26 s )\e)(e\ 

n'a'* 



a|e)(g| + n|e)(s| 



Mif 



H.c. 



(1) 



Here, the Rabi frequencies of the signal and the Stokes 
fields are a = £fj,/h and a' — S'fx'/h, where £ and £' are 
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the corresponding slowly-varying envelopes, and n and 
\J are the (real) dipole matrix elements of the respective 
transitions. 

In the undepleted pump and adiabatic approximations, 



the Fourier components of the signal field £(u) and the 
Stokes field £'(co) propagate according to the coupled dif- 
ferential equations, to linear order in a and a' [TBI [T51 123] : 



£(z,w) 
£'*{z,u) 



dry 



S-Ss + u + i^o -tt 2 /A hi 



n* 2 /A 



hf 



£(z,w) 
£'*{z,uj) 



(2) 



where we take into account the optical polarization decay 

in'i 2 

rate 7 and the ground state decay rate 70 ~ 7™ + 7 L at~ 

[25] . and have set F = \fl\ 2 + [7 - i(5 - 28 s + uj)][y - 
i(S — 5 S + lo)]. Under the assumption that CI is real 
and z-independent, Eq. ^ can be solved analytically for 
£{z,lo) and £'*(z, ui), for the conditions corresponding 
to our experiment: £'*(0,u) = —f£(0,ui): the Stokes 
seed has the same initial temporal lincshape as the sig- 



nal pulse, but with an opposite phase and with some 
amplitude scaling factor < / < 1. Defining the Ra- 
man detuning as Ar = — f2 2 /Ahf, and with f3(u>) = 
v/fT^TtrfTT^F+iAX a(cu) = i*L(S-5 a+w + 

«7o), and £(u>) = |^j;/3(w), we find the following analytic 
expressions for the Fourier components of the signal and 
Stokes fields IEJ EB]: 



£(z,oj) = £(0,ci;)e <<r < w > 
£'*(z,uj) = -f£(0,u)e ia ^ 



cosh[£Hz] + i I - ) sinh[£(w)z] 



cosh[eH,]-,(g|-I^)sinh K H„ 



(3) 
(4) 



Eqs. ^ and Q fully describe the propagation of the light 
fields through the atomic medium. Theoretically, the 
measured transmission spectra of the signal and Stokes 
fields are computed as and |£'(L)|, respectively 

using Eqs. ([3| and Q. We interpret equations ^ and 
Q by first applying a few simplifications, similar to 
Ref. [25], by first shifting the two-photon detuning by 
the light shift, defining 6 = S — 5 S (setting u = for 
continuous wave measurements), and considering large 
|?| > 2|A R |, but also \S\ < ft, and \S\ > 70. Un- 
der these assumptions, ~ iS, and for our parameters, 
2£L = 2iaL 



'WW ~ [^/(AF ' whcrc thc denom " 
inators of the two terms represent, respectively, the in- 
verse of the EIT group delay and the square of the width 
of the EIT transparency window. Under these approxi- 
mations, which hold well in our experiments, the signal 
and Stokes amplitudes after the cell are: 



\£(L)\ = £(0) 
\£'(L)\ = £(0) 



o 2 
A M <5 



o 2 



A hf <5 



(1 - e MaL ) ~ f 



(5) 
(6) 



These expressions allow us to interpret the transmis- 
sion spectra for both fields in terms of an interference 



between EIT and FWM effects, where the FWM effect 
is represented by the terms proportional to l/Ay [26] . 
At small optical depth, the effect of FWM on the signal 
field transmission is negligible, and we observe a typical 
symmetric EIT transmission peak (not shown). As opti- 
cal depth increases, the FWM term (oc l/Ajj) in Eq. ^ 
becomes more noticeable relative to the EIT term, which 
reduces with d due to thc narrowing of thc EIT window. 



Since the phase of 



in Eq. 



S/(v g /L), for 



S > 0, S — WKVg/L gives destructive interference (and 
hence dips in the spectrum) for even n, and constructive 
interference for odd n. For S < 0, the opposite case is 
true: even n yields constructive interference; odd n yields 
destructive interference. 



IV. SPECTRAL MEASUREMENTS 

We record the transmission spectra for continuous sig- 
nal and Stokes fields by sweeping the EOM frequency, 
which simultaneously changes the two-photon detuning 
8, and measuring their amplitude variations after the cell. 
Fig. [2ja,b) depicts the experimental results for an opti- 
cal depth of 2c? = 52. Solid lines represent the spectra 
corresponding to the signal field transmission; dashed 
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FIG. 2: (Color Online) (a, b) Signal and stokes amplitude spectra at an optical depth of 2d = 52 with a beam diameter of 
4 mm and £l/2ir = 9 MHz. Black traces are with a full Stokes seed present. Red traces are with the Stokes seed amplitude 
attenuated to ~ V0.05 of the signal field's amplitude, (a', b') Corresponding theoretical predictions, (c, d) Signal and stokes 
amplitude spectra at an optical depth of 2d = 98. (c', d') Corresponding theoretical predictions. 




FIG. 3: (Color Online) (a, b) Signal and stokes amplitude spectra at an optical depth of 2d — 52 with a beam diameter of 2.6 
mm and 0,/2-k = 14 MHz. Black traces are with a full Stokes seed present. Red traces are with the Stokes seed amplitude 
attenuated to ~ V0.05 of the signal field's amplitude, (a', b') Corresponding theoretical predictions, (c, d) Signal and stokes 
amplitude spectra at an optical depth of 2d = 110. (c', d') Corresponding theoretical predictions. 



lines correspond to the Stokes field transmission spec- 
tra. Black lines are with no Stokes filtering (i.e., / = 1); 
red lines are with a Stokes intensity attenuation so that 
/ = V0.05. These two values of / (1 and V0.05) are 
shown by horizontal dashed blue lines in Figs.[2jb,b') and 
represent the input Stokes amplitude. With the reduced 
Stokes seed (red curves), the effects of FWM are sup- 
pressed, and we observe a slightly amplified and nearly 
symmetric EIT transmission peak. However, when the 
full Stokes seed field is present (black curves), on one 
hand we observe more gain in the signal field, but on the 
other hand the FWM/EIT destructive interference be- 



comes more evident by the presence of a "knee" in the 
signal spectra for small negative detunings. These results 
are in very good agreement with the predictions of the 
theory [see Fig. j2)(a',b% which are calculated from the 
full expressions in Eqs. ^ and Q with no free param- 
eters, where f2, 70, d, and 7 were computed as in Ref. 

El- 

Spectra taken at higher optical depth reveal more clear 
evidence of the constructive and destructive interference 
between EIT and FWM. Fig.[2^c,d) presents similar spec- 
tra recorded for the same signal and Stokes fields, but 
at an optical depth of 2c? = 98. In the case of no 
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Stokes attenuation (black curves), the theoretical curve 
in Fig. [2^c') exhibits, as expected, destructive interfer- 
ence at 5 = 5 — 5 S = nirVg/L for n = —3,-1, and 2 
[here ttv 9 /L = (2tt)31 kHz and S s = (2tt)36 kHz]. While 
slightly shifted, these three points of destructive inter- 
ference are also clearly visible in the experimental mea- 
surement of Fig. |2jc). Even when the Stokes seed field 
is suppressed, its presence leads to significant distortions 
in the signal transmission resonance. 

We repeated similar spectral measurements after re- 
ducing the diameter of the beam by a factor of 1.5, which 
increased the control field Rabi frequency to H/2tt = 14 
MHz, and corresponded to a larger light shift of 5 s /2ir ss 
85 kHz. Larger control intensity and smaller beam size 
allowed us to reduce the control field absorption at high 
optical depths and stay within the theoretical model's 
approximations. Fig. [3] shows the experimental and cor- 
responding theoretical spectra for optical depths of 52 
and 110. The larger Rabi frequency results in a larger v g 
than above, and thus more closely spaced spectral dips 
and peaks. There is an excellent agreement between the 
experiment and the theory for the Stokes spectra, and 
for the signal spectra at negative detuning. However, the 
theoretical model for signal transmission diverges from 
experimental observations at positive detuning, indicat- 
ing the presence of some unaccounted mechanisms such 
as nonunity control field refractive index, atomic diffu- 
sion [37], and/or the effects of the multi-level structure 
of the atoms. 



V. SLOW LIGHT MEASUREMENTS 

In this section we discuss the slow light regime for the 
signal field pulses in the presence of the co-propagating 
seeded Stokes field. In particular, we are interested in 
the prospect of manipulating the signal pulse group de- 
lay and amplitude via the controllable amplitude of the 
input Stokes field. It is convenient to use Eqs. pi, Q to 
analyze the dynamics of each pulse propagation through 
the medium by calculating the variation in the com- 
plex amplitudes of both fields for each spectral compo- 
nent of the input pulse and then by Fourier transform- 
ing the resulting expressions back into the temporal do- 
main. The group delay of the signal field is determined 
from the acquired phase, which consists of two contribu- 
tions. The first from the first exponential in Eq.([3]), and 
it is the same for all spectral components of the pulse: 
t = j^~Re[a(u>)z) ss ^To 5 "- Notably, this value is ex- 
actly half of the pulse delay expected from the pure EIT 
system. The second contribution is from the expression 
in brackets in Eq. The value of this additional de- 
lay depends explicitly on the detuning of the signal pulse 
from resonance and may vary significantly for different 
spectral components of the pulse. Below we discussed 
three distinct scenarios for the pulse two-photon detun- 
ing: 5 = 2|A R |, 5 = 25 s , and 5 = 0. 

Figs. |4ja,b) and |5ja,b) correspondingly present the 



experimental data for 6 /xs-long (FWHM) signal and 
Stokes pulses (which corresponds to a bandwidth of 
±(27r)31 kHz around the carrier frequency) when the full 
Stokes field is present, and when the Stokes field is sup- 
pressed. Respective graphs (a') and (b') give the predic- 
tion of the calculations based on the complete solutions of 
Eqs. ([3]),(|4]). In these calculations, we use a control field 
with Rabi frequency £l/2w = 14 MHz, corresponding to 
5 s /2n = 84 kHz and A r /2tt = -28 kHz. 

For more insight into the spectral dynamics of the 
pulse, we also plot the calculated time delay experienced 
by the signal field spectral components u [Figs. |4jc) and 
gc)], and the spectral gain \£(w, L)\/\£(uj, 0)| [Figs.gd) 
and [5^d)]. These last graphs also show the spectral 
bandwidth of the input pulse (the blue, dashed curve) 
for reference. 

Case I: 6= 2|A R | 

The data corresponding to the case 5 = 2|Ar| = 
(27r)56 kHz are shown in Figs. [4 and [5] in red. Although 
the theoretical calculations predict only a weak depen- 
dence of the signal pulse on the amplitude of the Stokes 
seed, the experimental data show much stronger depen- 
dence: for the full seeded Stokes field the experimen- 
tal pulse shows small delay and noticeable attenuation, 
while when the Stokes field is small, the signal pulse ex- 
periences some small gain and much larger delay. The 
latter is much closer to the theoretical expectations of 
gain of ss 1.5 and the delay of ss 6/xs. The experimental 
results for the Stokes field, however, match the theory 
much more accurately, which may indicate that the ab- 
sorption of the signal field is undcrrepresented by the 
model. 

The analysis of the spectral gain and delay for both 
cases (/ = 1 and / = y/0.05) provides some qualita- 
tive understanding of the observed pulse behavior. For 
the case of the full Stokes field, Fig. Qc) shows that all 
the signal spectral components with u> > experience a 
roughly uniform delay of ss 5^ts, whereas the components 
lu < experience a longer delay, resulting in the pulse 
spreading as it traverses the cell. However, this effect is 
somewhat suppressed by lower spectral gain for lu < 
(Fig. [IJi). With the Stokes field attenuated [Fig. ^c,d)), 
we expect that all signal spectral components lu experi- 
ence a uniform delay of ss 6fis, and a slight gain, resulting 
in delayed propagation with little pulscshape distortion, 
as corroborated well in Fig. (5ja). 

Case II: 5 = 2S S 

The green curves in Figs. [4] and [5] depict the results 
of slow light experiments with a two-photon detuning of 
5 = 25 s = (2tt)168 kHz. Figs. Rc,d) illustrate that for 
the unfiltered Stokes field (/ = 1), all frequency compo- 
nents of the initial signal pulse experience a nearly iden- 
tical delay of ss 5[is — indicating very little pulse spread. 
Simultaneously, the central component should be ampli- 
fied by a factor of ss 1.8. This prediction matches well 
with both the experimental [Figs. [3Ja,b)] and theoretical 
[Figs. Qa'jb')] pulses. When the Stokes seed is attenu- 
ated (/ = \/0.05), the signal pulse experiences a slightly 
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FIG. 4: (Color Online) (a) Slow light on the signal channel at an optical depth of 2d = 110 with no Stokes seed attenuation 
(i.e., f — 1) for signal detunings 8 = 2|Ar| (red curves), S = 2<5 3 (green curves), and 5 = (black curves). The thin blue curve 
is the initial reference pulse, (b) Corresponding Stokes channel, (a', b') Corresponding theoretical predictions from Eqs {3},|4). 
(c) Theoretical total delay dispersion experienced by the signal pulse frequency components after traversing length L. The blue 
dashed curve depicts the frequency spread of the input signal pulse, (d) Predicted signal gain dispersion. 
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FIG. 5: (Color Online) Same as Fig. [4] but with Stokes seed attenuation (/ 



longer delay of ss 5.8/is, but will also be less amplified, 
according to Fig. [5Fd). The experimental result [green 
curve in Fig. [5j a)] reproduces this predicted delay, but 
shows a small attenuation rather than gain, possibly in- 
dicating the presence of an additional decay mechanism. 

Case III: 6 = The black curves in Figs. [4] and [5] 
depict the results of slow light experiments with a two- 
photon detuning of 6 = 0. This case most clearly demon- 
strates the merits of Stokes seed attenuation. For the 
unfiltered Stokes seed / = 1, different spectral compo- 



nents will acquire very different phase and gain while 
propagating through the interaction region. In particu- 
lar, Fig |4^c) shows large variation in the spectral delay 
- from 14/xs for central frequencies to 5^s for the far- 
ther detuned components. When combined with the gain 
curve, shown in Fig. |3Jd), such variations should greatly 
distort the shape of the output pulses. In fact, the ex- 
pected output closely resembles a double-peaked pulse, 
and is quite similar to that observed in the experiment. 
When the Stokes seed is filtered (/ = V0.05), as in Fig. [H] 
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the differential delay is suppressed, and all spectral com- 
ponents experience a common delay of nearly 7/zs, but 
at the sacrifice of gain, which is < 1. Figs. [5][a,b) show 
excellent agreement with the corresponding theoretical 
curves. 

VI. CONCLUSIONS 

In conclusion, we have demonstrated that both steady- 
state and dynamic properties of the signal field propa- 
gating under the EIT conditions are strongly effected by 
resonant four-wave mixing that arise under the condi- 
tions of EIT at high optical depth. This process is well- 



modeled by a simple double-A system, where the output 
signal and stokes field amplitudes are the results of inter- 
ference of "traditional" EIT and FWM. We have shown 
that by attenuating the amplitude of the seeded Stokes 
field, we can partially control the optical properties of 
the medium for the signal field. Moreover, by adjust- 
ing the central frequency of the input signal field around 
two-photon resonance, in the presence of Stokes seed field 
we can achieve longer pulse delay and/or amplification of 
the signal pulse. 

The authors thank M. D. Lukin for useful discus- 
sions. This research was supported by NSF grant PHY- 
0758010, Jeffress Research grant J-847, and the College 
of William & Mary. 
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